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Abstract 

We derive semiclassical periodic orbit expansions for a correlation function of the Wigner 
time delay. We consider the Fourier transform of the two-point correlation function, the form 
factor K(t, x, y, M), that depends on the number of open channels M, a non-symmetry breaking 
parameter x, and a symmetry breaking parameter y. Several terms in the Taylor expansion 
about r = 0, which depend on all parameters, are shown to be identical to those obtained from 
Random Matrix Theory. 
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1 Introduction 



Random Matrix Theory (RMT) and semiclassical approximations are two alternative theoret- 
ical approaches for quantum systems that are chaotic in their classical limit. RMT describes 
universal, system-independent, statistical quantum fluctuations by modelling the Hamiltonian 
with matrices that have random elements [1, 2, 3]. Semiclassical methods on the other hand 
yield approximations in terms of classical trajectories of individual chaotic systems, asymptot- 
ically valid as h — > 0, and can also describe non-universal system-specific quantum properties 
[4, 5]. In recent years, however, a main emphasis in the application of semiclassical methods has 
been to obtain the universal results of RMT from the asymptotic properties of long trajectories 
and thus provide a theoretical foundation for the applicability of RMT in chaotic systems. 

The first quantities that were treated in this way were spectral correlation functions of 
closed systems which have a semiclassical expansion in terms of a double sum over periodic 
orbits. For the spectral form factor K(r), the Fourier transform of the two-point correlation 
function of the density of states, the terms in the Taylor expansion about r = were derived 
semiclassically in agreement with RMT. The first term in the Taylor expansion was obtained 
from the diagonal approximation for the double sum over periodic orbits which pairs orbits 
with themselves or their time reverse [6, 7]. Higher order terms are due to pairs of different, 
but strongly correlated periodic orbits. The first off-diagonal term, i.e. the second term in the 
Taylor expansion, was calculated for uniformly hyperbolic systems in [8, 9]. The third order for 
general chaotic systems was obtained in [10], and all higher order terms in [11, 12, 13]. Similar 
methods have been applied to derive off-diagonal terms, for example, for the conductance 
[14, 15], the GOE-GUE transition [16], the shot-noise [17] and parametric correlations [18, 19]. 

In this article we consider two-point correlation functions of the Wigner time delay in 
open systems, for which there are very comprehensive results from RMT available [20, 21, 
22]. We evaluate a semiclassical periodic orbit expansion for the Fourier transform of the 
two-point correlation function, the form factor K(r,x,y, M). It depends on a (scaled) non- 
symmetry breaking parameter difference x, the number of open channels M, and parameter 
y that depends on a time reversal symmetry breaking magnetic field. Earlier work treated 
K(t, M) and K(r, y, M) in the diagonal approximation [23, 24]. We derive several terms in the 
Taylor expansion of K(t, x, y, M) around r = and compare them to RMT. We find complete 
agreement for the terms in this expansion that depend on the four parameters r, x, y, and M. 

In section 2 we introduce the form factor K(r,x,y, M), and in section 3 we state results 
of RMT for it. Sections 4 and 5 treat the case without symmetry breaking parameter y and 
contain derivations of diagonal and off-diagonal terms for K(t, x, M), respectively. The general 
case with additional parameter y is treated in section 6, and section 7 contains our conclusions. 



2 The form factor for the time delay 



The Wigner time delay is a measure of the typical time that is spent in a scattering process in 
comparison to free motion. It is defined in terms of the scattering matrix S(E) as [25, 26] 



r W (E) = -§TV 



S*(E)—S(E) 



MdE v ; ' 



(1) 



where M is the number of open scattering channels. We allow the system to depend on external 
parameters X and Y and will include this later in the notation. A comprehensive overview of 
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the history and the applications of the time delay is given in the review [27]. 

One important property of the Wigner time delay is that it can be related to a density of 
states [28, 27] 

rw(E) = 2 -^d(E) « ^ (d(E) + d™(E)) (2) 

The density d(E) here is the difference between the level density of the scattering system and a 
free or reference system, because the level density itself is infinite. In the semiclassical regime 
it is divided into a mean part d(E), and an oscillatory part d osc (E). The first term in (2) 
represents the average time delay which can be expressed as T H /M where T H = 2irhd(E) is the 
Heisenberg time. We will discuss later that the oscillatory part is semiclassically approximated 
by the periodic orbits in the repeller. In the following we will consider the case of a two- 
dimensional open chaotic cavity with attached leads. Then d(E) is given by Weyl's law for the 
closed cavity d(E) ~ Q/(2nh) 2 , where Q is the volume of the surface of constant energy E in 
phase space, and d osc (E) is given in terms of the periodic orbits in the open cavity [24]. 

Quantum fluctuations of the Wigner time delay are characterized by its correlation functions. 
If one compares the time delay at different energy and parameter values one has to scale the 
energy difference as well as the parameter difference in order to obtain universal correlation 
functions in chaotic systems. Similarly as for closed systems one can define a new energy 
parameter E by 

E = N{E,X), (3) 

where N(E,X) is the mean part of the integrated level density N(E,X) = J_ dE' d(E',X). 
In term of this 'unfolded' energy variable the level density has a mean value of one. Small 
intervals of the two energies are related by AE ss AE/d(E). In closed systems, the external 
parameter X is unfolded by requiring that the variance of the level velocities dE n /dX, of 
the energy levels E n , is unity in terms of the new parameters X and E. One then obtains 
AX ~ AX /a, where a 2 is the variance of the level velocities. For open systems we scale the 
parameter difference again by a quantity a whose value we will discuss later (see equation (29)). 

In the semiclassical limit one can then define a universal correlation function by (see [19] 
for closed systems) 

, ,\ (^(E + % + %,x + £,Y) ^(jg-a-fg.x-^y) ) 

R 2 {r],x,y,M) _ - (4) 

where y is a function of Y to be specified later. The average is carried out over an energy 
interval AE, satisfying E <C AE <C l/d(E), and over a parameter interval AX. We note that 
changes of order one of x and rj correspond to very small changes on the classical scale. The 
term xp/2a accounts for the change of the energy when X is changed while keeping E fixed 



BE 
dX 



ON/OX (5) 



& dN/dE 

We can also express the two-point correlation function in terms of the level density 

„. , ( d ° SC (E+^+^X + i) d osc (E - £ - g ,X - £) ) 

RM 1(E-Xf (6) 
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The parametric form factor for the time delay is obtained by a Fourier transform of this 
two-point correlation function 

/oo 
R 2 (r],x)e- 2mr!T dr) . (7) 
-oo 

3 Results from Random Matrix Theory 

In this section we review the RMT calculations for the form factor for the time delay. We are 
interested in the small r expansion of the form factor and we restrict ourselves to the case r < 1 
in the following. We consider first the results for the case where we do not have a parameter that 
breaks time reversal symmetry, and state the results for the Gaussian Orthogonal Ensemble 
(GOE) and the Gaussian Unitary Ensemble (GUE) that are relevant for systems with and 
without time reversal symmetry, respectively. The semiclassical result will be derived for perfect 
transmission T = 1 in all scattering channels, and we state the RMT results only for this case. 

The parametric two-point correlation function i? 2 (^7, M) for the GUE case was derived in 
[21] and is given by 

R^ E ( V , x,M)= 1 - £ dA dAx cosMA! - A)) (i±J* ) " e -*>W-#)/2 (8) 

After performing the Fourier transform in (7) to obtain the form factor we arrive at 

if GUE (r, x, M) = i J 1 dA dAi <f(Ai - A - 2r) (j^j " ^V*'-*)/* ( 9 ) 

In this equation we have removed a second delta function which does not contribute due to the 
restrictions r > and Ai > A. The remaining delta function provides the relation 2r = Ai — A. 
In the case r < 1, the domain of integration for Ai is reduced to 1 < Ai < 1 + 2r, and we obtain 

K GUE (r, x, M) = - [ 1+2T dA x ( 1 + Al ~ 2r ) V^) , r < 1 . (10) 
2 Ji V 1 + Ai / 

To get a series expansion of this integral, it is useful to remove the r dependence from the 
limits by a change of variable. We do so by using Ai = 1 + ryi, and we also define B = 2n 2 x 2 / k 
where k — 1 and 2 for the GUE and GOE cases, respectively. The integral then becomes 

iT GUE (r, x, M) = -e~ Br / dyi ( 1 - e^ 2 ^ (11) 

2 Jo V 2 + ryi/ 

Now we can expand the integrand in a power series in r. When we perform the integration, 
term by term, and extract an exponential we get the following expansion up to 9th order (by 
using MAPLE) 



K GVE (r,x,M)=e^ B+M ^ 



M 4 (2B-M) 2 5 M 6 (BM M 2 \ 7 

— t H r t — T 

6 24 15 V 15 20 



(1M(2B - Mf M\ 
V 720 + 28/ 



r ' 



(2B — M) 4 BM 401M 2 \ 
1921) 28~ + 10080 



(12) 
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The form factor for cases without parametric variation is obtained by setting B = (or equiv- 
alently x — 0). It is therefore 



,,,,, . . Mr 4 M 2 r 5 Mr 6 MV 7M 3 r 8 

^GUE M = e -Mr + + 

6 24 15 20 720 

Mr 8 M 4 r 9 401M 2 r 9 



+ "TTTTTTT + + 



28 1920 10080 



(13) 



For the GOE case, the correlation function can be found for the case without parametric 
variations in [20] , and for the case with parametric variations in [22] The parametric correlation 
function, which we consider here, is expressed by a triple integral 

R GOE ( V ,x,M) = TdA f °dX 1 f°dA 2 cos(tt7 ? (A-A 1 A 2 ))- (1 " A 2 )(A - A x A 2 ) 2 



(2AA 1 A 2 -A 2 -A 1 2 -A 2 2 + 1) 2 

x ( 1 + A \ M c -A 2 (2A 1 2 A 2 2 -A 2 -A 1 2 -A 2 2 +D/4 ^ 

\ Ai + A 2 y 

To obtain the parametric form factor we take the Fourier transform and obtain a sum of two 
delta functions. Again, because r is positive and AiA 2 > A only one of the delta function 
contributes, and we only include this one in the expression 

K GOE (r,x,M) = f dA TdAi /°°dA 2 5(A-A 1 A 2 -2r) (1 ~ A ' )(A ~ 

1 ' J-i Ji Ji 1 (2AAiA 2 - A 2 - Ai - A 2 + 1) 2 



1 + A 



M 



x / t 1 " j e -7r 2 a; 2 (2Ai 2 A2 2 -A 2 -Ai 2 -A 2 "+l)/4 

\ Ai + A 2 / 

The delta function yields the relation A = AiA 2 — 2r. As r < 1, our domain of integration for 
the other two variables is given by 1 < Ai < 1 + 2r and 1 < A 2 < When the integral over 

A is performed, we are left with 

K {r,x,U) = l A* I d ^ (1 + VV _ V _ V _ 4TT ( Al + A2 J 

x e -7r 2 x 2 (l+A 1 2 A 2 2 -A 1 2 -A 2 2 +4rA 1 A 2 -4T 2 )/4 ^g-j 

In order to compare to semiclassical results this integral is evaluated as a series in r. To do 
so we first remove the r dependence from the limits by introducing new variables defined by 
Ai = 1 + ryi and AiA 2 = 1 + ry 2 . After this change of variables, the integrand is expanded as a 
series in r. We also replace 7r 2 a; 2 by B (as defined after equation (11)), though we only include 
the first two orders for clarity 



K GOE (r,x,M) = fd yi f 
Jo Jyi 



2(1 - ym + yi 



yi (y 2 - 2) (4 - Vl y 2 + 2y x 2 - y 2 2 ) 
2{l-y 1 y 2 + y l 2 f 

r 2 + ...l (17) 



(2-y 2 )(2 + 2B + 2M -Q yi -y 2 ) 

4(1- ym + yi 2 ) 2 

Once the integrals are evaluated we get a series expansion for the parametric form factor for 
the time delay. In order to compare with semiclassical results, we again extract an exponential 
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term, to get the following result to seventh order 



K GOE (r,x,M)=e^ B+M ^ 



+ 



It - 2r 2 - (2B - M - 2)r 3 + ( 2B 
5(25 - M) 2 8B ISM 



7M 8 



3 IT 



12 



3 + 3 



+ 4 r 



55 2 WBM 17M 2 
~3 3 + 12 

41(25 -M) 3 115 2 
360 5~ 



45 + 



39M 32 



+ 7BM- 



5 

43M 2 
12 



+ 



325 212M 32 



t' + 



5 15 3 

By setting B = we get the form factor for the case without parametric correlations 



(18) 



K 



GOE, 



T,M) 



-Mr 



2r - 2r 2 + (M + 2)r 3 - 



+ 
+ 



5M 2 13M 

1 h 4 r 5 

12 3 



3 3y 

/17M 2 39M 32 



V 12 

/41M 3 43M 2 . 212M . 32 
V 360 + ~Y1 



+ 



+ 



15 



t' +... 



(19) 



Finally, we include the dependence on a further parameter y that breaks time-reversal 
symmetry and thus leads to a GOE-GUE transition as y is varied from zero to infinity. The 
RMT result can again be found in reference [22], the equation (2) in this article includes even 
more general cases. For our case, the correlation function has the form of the integral (14) with 
an additional factor in the integrand that is given by 



G(y) = exp (|(A 2 + 1-2A 2 )) 



x 



where 



[I - A 2 )(l + ^) cosh a + (A 2 - A 2 ) sinha + ^(2A 2 + A 2 - 1) sinha 



R = A 2 + A 2 . + A 2 - 2AAiA 2 - 1 , and a = -(1 - A 2 ) 



(20) 



(21) 



and we choose the parameters y\ and y 2 in [22] as y\ = y\ = y/8. The steps to get the form 
factor are exactly the same as for the GOE case and we give only the final result 

K(r,x,y,M) =e^ B+M ^ 

1 



2 T -{2 + y)r 2 +(^ + 2y+ l -y 2 -2B + M^ r 3 - ^ + 2y 



+V + 7V — 2B — By + -M + -My r 4 + 4 + 2y + -y l + -y d + —y 



6 



6 1 



24' 



~B + yM + ^(25 - M) 2 - ±(B - M)y - ^(2B - M)y 2 ) r 5 + . . . 



(22) 
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Unfortunately, one cannot see the transition from GOE to GUE in this formula. As will be- 
come clear from the semiclassical calculation, there are exponential terms of the form exp(— yr) 
which describe this transition, but these have been expanded when we expanded the integrand 
in powers of r. We tried to replace y by y/r in order to extract these exponential terms, but 
then we could not bring MAPLE to perform the integrals. However the equation (22) can be 
compared to the semiclassical result if we expand the exponentials there as well. 



4 Semiclassical approximation 

In this section we derive a semiclassical expression for the form factor. The oscillatory part of 
the density in (2) has exactly the same semiclassical form as in Gutzwiller's trace formula for 
the density of states [29, 30] 

d osc (E,X) « ^-Re VA^exp ( ^S 7 ] , where A, = ^ . (23) 

V ^ / /2 7V / |det(M 7 -l)| 

The sum runs over all periodic orbits in the repeller with period T 7 , repetition number i? 7 , 
stability matrix M 7 and Maslov index // 7 . 

In order to evaluate the two-point correlation function of the form factor we expand the 
action in first order in the energy difference and the parameter difference 

( E ± a ± x ± B K w x) ± T ^ E - x) Yi ± Q ^ E - x) To (24) 



where 



dX 



, = "a¥ + af (25) 



is the parametric velocity. In this and the next chapter we consider systems without a parameter 
that breaks time reversal symmetry. The case with such a parameter will be discussed in 
section 6. 

After inserting (23) and (24) into (4) and (7) and evaluating the integral to leading semi- 
classical order, we arrive at 

1 I T-^ HS 7 -S ,) \x(Q 1+ Q ,) / 71 + T,v/\\ 

K(r, x) = jr(Y, A J A;,e^^e^^5 [t - ±±±) J (26) 

where T H = 2nhd(E) is the Heisenberg time, and r = T/T H . Terms which have a sum of 
the actions in the exponent have been neglected, because they are averaged away. Equation 
(26) has exactly the same form as for the parametric form factor for closed systems [19]. The 
difference is that the sum is over the periodic orbits in an open system which obey a different 
sum rule than the periodic orbits in a closed system. 

The diagonal approximation involves pairs of orbits that are either identical or related by 
time reversal. It has the form 

K(r,x,M) = ^ (^£\A^e^5(T - , (27) 

where k is 2 if the system has time reversal symmetry and 1 if it does not. 
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An important input in the following semiclassical calculation is the distribution of the 
parametric velocities Q 7 in the limit of very long periodic orbits. The Q 7 have a mean value of 
zero and a variance proportional to their period [31] 

(Q 7 ) = , (Q 7 > ~ aT , T -> oo (28) 

where the averages are performed over trajectories with period around T. 

The main assumption that we will use in the following is that the Q 7 have a Gaussian 
distribution for long periodic orbits. This has been motivated e.g. in [32]. In closed systems 
the proportionality factor a in (28) is related to the variance of the level velocities by [33, 34] 

- 2 ~fi- < 29 > 

We will see that we obtain the universal form factor when we use this relation for the scaling 
factor a in (4) in open systems as well. 

With the Gaussian assumption for the distribution of the parametric velocities and (29) the 
average over the Q 7 can be performed, assuming that it can be done independently, and we 
obtain 

, ix Ql . X 2 aT D rp irp 

( e — ) = e"^^ = q-bt/Th ^ 

where B = 2n 2 x 2 / 'k has been introduced after equation (11). The remaining sum over periodic 
orbits can be evaluated with a modification of the Hannay-Ozorio de Almeida sum rule [6]. 
This modification takes into account the probability that classical trajectories escape in open 
systems. The classical escape rate, or the inverse of the mean time spent in the cavity, can be 
expressed as fj, — M/Th (see e.g. [35]), and the sum rule takes the form [36] 

£ |A 7 | 2 5(T - T 7 ) » Te- MT ' TH . (31) 

7 

We find that the diagonal approximation is given by 

K d (r,x,M) = ktq- { - b+m)t (32) 

in agreement with the first term in the expansion of the random matrix results, (12) and 
(18). Previous evaluations of the diagonal approximation of the form factor for the time delay, 
respectively its Fourier transform, were performed in [23, 24]. The second reference also treats 
time reversal symmetry breaking. 



5 Off-diagonal contributions 

The correlations between different periodic orbits in open systems that are responsible for the 
off-diagonal terms in the expansion of the form factor have the same origin as in closed systems 
[8, 11, 12]. The general formalism for evaluating the correlations between orbits for arbitrary 
order has been developed in [11, 12, 13]. 

Long periodic orbits have many close self-encounters in which two or more stretches of an 
orbit are almost identical, possibly up to time reversal. These orbit stretches are connected by 
long parts of the orbit which are called "loops". There exist other periodic orbits which are 
almost identical along the loops, but they differ in the way in which the loops are connected 
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in the encounter regions. As a consequence the actions, periods and stabilities of these orbits 
are strongly correlated. One can classify correlated orbit pairs according to "structures" that 
are characterized by the number of encounter regions V in which the loops are connected in 
a different way, the number of involved orbit stretches l a in each encounter region a, and the 
way in which the loops are connected by these stretches. Structures can be put in a one-to-one 
relation with permutation matrices that describe the reconnections of the loops in the encounter 
region. It is a combinatorial problem to find all possible structures that describe correlations 
between periodic orbits. 

The number and types of encounter regions are specified by a vector v whose Z-th component, 
vi, denotes the number of encounter regions with I stretches. The total number of orbit stretches 
is L. Hence 

l>2 a l>2 

In general there are many structures with the same vector v, because, for example, there are 
many different ways in which the loops can be reconnected in the encounter regions, and the 
number of structures with the same vector v is denoted by Af(v). 

The summation over all orbit pairs with the same structure can be performed by using 
that long periodic orbits are uniformly distributed over the surface of constant energy in closed 
systems. In open systems one has take additionally into account that trajectories can escape 
to infinity. One combines the contributions from all structures with the same vector v, and 
the semiclassical form factor is then obtained by summing over all v and adding the diagonal 
approximation 

K sc (r) = kt + Y^Kv(t) (34) 

V 

where K v (t) is the contribution to the form factor from all structures with the same v. 

The derivation of K v (t) for the spectral form factor can be found in [12, 13]. We state here 
only those details that are relevant for the following calculations. The semiclassical expression 
for K v (t) can be expressed as 

fixed v p / \ 

Kv (t) = —Y \A,\ 2 e iAS ^ h 5(T - T 7 ) = N(v)kt / d L ~ v s d L ~ v u Wt ^ U > j™' h . (35) 
Th (7,70 J L 

Here the components of the vectors s and u are coordinates along the stable and unstable 
manifolds in Poincare sections in the V encounter regions, that describe the relative positions 
of the orbit stretches. In (35) the amplitudes and periods of the two correlated orbits are set 
equal. wt(s, u) is the density of self-encounters for a given structure and separation coordinates 
s and u. For long orbits it is given asymptotically by 

L L\W-vi[ Ac ■ 1 ) 

Here a labels the V different encounters, each being a /^-encounter, and t^ nc is the time for the 
traversal of the encounter region a. The traversal time is a function of s and u. 

The only property that is needed for an evaluation of the integral in (35) is the semiclassical 
relation 



J n dsa J 




du aj (C)V^ ^ - { 11 - ^ ° r k - 1 . (37) 
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This considerably simplifies the evaluation of the contribution K v (t), because after an expan- 
sion of the numerator of (36) the only terms that survive are the ones that contain a product 
of all encounter times t^ nc which is cancelled by the denominator. In this way the contribution 
K v {t) was evaluated for the spectral form factor in [11, 12]. By deriving a recursion relation 
involving the numbers J\f(v) the semiclassical expansion could be summed up, and it was shown 
that it agrees with form factor of RMT for r < 1. 



5.1 Off-diagonal terms without parametric correlations 

Let us now come back to the form factor for the time delay. We consider first the case without 
parametric correlations, i.e. x = 0. We have to evaluate 

K v (t, M) = —J2 \ A i\ e s 6 ( T - T i) ( 38 ) 

(7:7') 

where the sum is over all correlated pairs of periodic orbits in open systems with the same 
number and type of encounter regions, specified by v. The density wt(s, u) in (35) has to be 
modified by the survival probability of long classical trajectories in open systems, similarly as 
was done in the sum rule for the diagonal approximation (31). If there would be no systematic 
correlations between different parts of an orbit then this survival property would be exp(— fiT). 
However the density wt{s,u) describes orbits that all have the same number and type of 
encounter regions. In an encounter region there are different stretches of an orbit which are 
almost identical. If an orbit does not escape during the traversal of one stretch in an encounter 
region then it has a negligible probability of escaping during the traversal of the other stretches 
in the same encounter region. Hence for each encounter region the survival time T should be 
reduced by (l a — l)t" nc . The total survival probability is then 



(39) 



This is analogous to the survival probability for correlated pairs of open trajectories that con- 
tribute to the conductance [15]. 

To summarize, the contribution of all orbit pairs with encounter regions described by v is 

fixed v r / \ 

K vM {t) = ±- Y \A,\ 2 e iAs ^ h 5(T H r -T 7 ) = N{v)kt / d L ~ v s d L ~ v u Wt[S > U) e isu / h (40) 
Th J L 

(7.7 ) 

where the density of self-encounters is now 

W T (S,U) _ e-^/^T^-EJaCj^'na^- 1 ^^ 



(41) 



which has been obtained by multiplying (36) by (39). Again because of (37) the only terms 
that contribute in the semiclassical limit are those where the encounter times in the numerator 
and denominator cancel exactly. As a first step we can expand the exponentials as a power 
series up to first order 

wtM . e- MT ' T »T(T - E Q 11 Jl + (L ~ l)Mt a enc /T H ) . , 
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The method to obtain K v (r, M) from the expression (42) is straighforward. For a given vector 
v (which specifies the set of the l a ) the numerator is expanded in all encounter times t" nc and 
only those terms are kept for which the product of the encounter times is exactly cancelled 
by that in the denominator. This is conveniently done with a computer program. K v (t, M) 
follows by multiplying the result with J\f(v)K,r(27rh) L ~ v , because of (40) and (37). Note that in 
order to obtain a product of the V different encounter times in the numerator we can take r of 
them from the product over a and V — r of them from the bracket with the power L — l, where 
< r < V. This gives a term to the form factor which is proportional to q- Mt t l ~ v+1+t M r . 

Table 1 shows the result for the different types of encounters with L — V < 4. The vectors 
v are represented in the form (2) V2 (3)'° 3 . . . and the horizontal lines separate vectors v with 
different value of L — V. The numbers N(v) can be calculated by combinatorial methods [13]. 
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Table 1: Contribution of different types of orbit pairs to the form factor for the time delay 



To find the total contribution to the form factor we now multiply the middle column that 
contains K v (t,M)/(kN(v)) by k and N(v), add the diagonal approximation and sum over 
different v. If we do that for all orbits pairs with L — V < 8 for the case without time reversal 
symmetry (k — 1), we obtain the expansion for the form factor in r up to 9th order 



K(t, M) = e 



-Mt 



Mt 4 M 2 t 5 
r — + 



2^7 



3^8 



7M t 
720~ 



24 
Mr 8 
^8~ 



Mi" Mr 
+ 



15 



M 4 r 9 
+ 1920 + 



20 

401MV 
10080 



+ 



(43) 



This agrees with the terms of the expansion (13) in the section on RMT. 

For systems with time reversal symmetry (k = 2) we sum over all contributions with L— V < 



11 



6, and obtain the expansion of the open form factor in r up to 7th order. 

'7M 8^ 



K(t, M) = e 



-Mr 



, 5M 2 13M \ - /17M 2 39M 32 \ 

+ 1 h 4 r 5 1 1 r 6 

'12 3 J V 12 5 5 



(44) 



,),-( 

/41M 3 43M 2 212M 32 \ 7 

+ {^w + — + — + T) T + 

This agrees with the expansion (19) of the GOE result. 

5.2 Off-diagonal terms with parametric correlations 

Let us now consider the parametric form factor K(t,x, M). We have to evaluate 

.. fixed v . , „ „ . 

K v (t,x,M) = — V lA/e^^e— 5 (T - T 7 ) (45) 

(7,7') 

It now contains an additional term involving the parametric velocities Q 7 . As for the diagonal 
approximation we assume that the average over the Q 7 can be performed independently from 
the actions and amplitudes of the orbits. However, we have to look again at the systematic 
correlations between different parts of the same periodic orbit that are implied by fixing the 
vector v [19]. In the encounter regions the different stretches of an orbit are almost identical, and 
the same holds for the changes of the action along these stretches when an external parameter 
is varied. Hence we should consider the average over the parametric velocities for the loops and 
encounter regions separately. The contribution from a loop is 



c 



*- ) = e -BT io o P /T H (46) 



while the contribution from the / orbit stretches in an /-encounter region is 

= e -Bl 2 t enc /T H _ ^ 



This means that the average over the parametric velocities is now given by 



e ^ \ = e -B(T~E a ^ nc )/T He -BZ a ll^ nc /T H ( 4 g) 



The contribution of all orbit pairs with encounter regions described by v then have the form 

fixed v „ / \ 

L-V„jL-V u ZT{S,U) eisu/h (4g) 



K v (r,x,M) = — V |A 7 |V A5 - /?l (5(T-T 7 ) = N{v)kt [ d L ~ v s d 



(7,7') 

were zt(s,u) is the product of the density of self-encounters wt(s,u) in open systems in (41) 
and the average over the phase factors in (48) 

z T (s,u) w T (s,u) i i-OiV e-( B+ ^ T / r *T(r-EaU^) L - 1 nae- (B,a -^ )(,a - 1) ^ /TH 

(50) 
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As before, because of (37), the only terms that contribute in the semiclassical limit are those 
where the encounter times in the numerator and denominator cancel exactly. As a first step 
we can expand the exponentials as a power series up to first order 

We proceed as in the previous subsection. For a given vector v the numerator is expanded in all 
encounter times t" nc and only those terms are kept for which the product of the encounter times 
is exactly cancelled by that in the denominator. This is done with a computer program. We 
multiply the result with N (v) kt(2tiK) l ~ v and obtain K v (t, x, M). Note that in order to obtain 
a product of the V different encounter times in the numerator we can take r of them from the 
product over a and V — r of them from the bracket with the power L — 1, where < r < V . 
This gives terms to the form factor which are proportional to e -( B + M ) T T L - v + 1 + r B s M r ~ s , where 
< s < r. 

The contribution of orbits for the different types of encounters are shown for orbit pairs 
with L — V < 4 in table 2. The double horizontal lines separate vectors v with different value 
of L - V. 
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Table 2: Contribution of different types of orbit pairs to the parametric form factor for the 
time delay 

To obtain the total contribution to the form factor we multiply the middle column that 
contains K v (t,x,M)/(kN(v)) by k and N(v), add the diagonal approximation and sum over 
different v. If we do that for all orbits pairs with L — V < 8 for the case without time reversal 
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symmetry (k — 1), we obtain the expansion for the form factor in r up to 9th order 



K{t,x,M) =e^ B+M ^ T 
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(52) 



This agrees with the expansion (12) in the chapter on RMT. If we remove the parametric 
correlations by setting B = we obtain the result (43). If we close the system by setting 
M = we get the expansion of the parametric form factor [19]. 

For systems with time reversal symmetry (/c = 2) we sum over all contributions with L — V < 
6, and obtain the expansion of the form factor in r up to 7th order. 



K(t,x, M) = e^ B+M ^ T 
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(53) 



This agrees with the expansion of the RMT result (18). For 5 = it agrees with the result 
(44) in the previous section, and for M = with the parametric form factor [19]. 



6 GOE-GUE transition 

We introduce now second parameter, a magnetic field Y, that takes a time reversal invariant 
system (at Y — 0) and breaks the symmetry as it is increased (Y — > oo). The semiclassical 
evaluation of the form factor for closed systems with a time-reversal symmetry breaking mag- 
netic field has been carried out in [37, 16, 18], and this calculation can be transferred to the 
parametric form factor for the time delay. 

The magnetic field adds a term 1 {B) = J'^Adq to the action, where A is the vector 
potential of the magnetic field that is assumed to be constant and perpendicular to our two- 
dimensional system. The semiclassical form factor is 

1 /_-^ i(S 7 -S ,) i(9 7 -9 ,) iz(Q 7 + Q ,) / T +T/\\ 

K(t, x, y, M) = — [t - -JLLjL J j (54) 

The sum is over all periodic orbits of the open system, and we have included the term with the 
parametric velocities that arise from a variation of the first parameter X. In the semiclassical 
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regime a very small magnetic field leads to a transition from from GOE to GUE statistics, and 
the influence of the magnetic field on the classical motion can be neglected. 

One ingredient in the following calculation is the distribution of the phases 9 7 (B). As for 
the parametric velocities it is assumed that they have a Gaussian distribution with zero mean 
and a variance proportional to the period of the orbit 

(0 7 ) = , (02) ~ 2Y 2 DT , T^oo (55) 

where the averages are performed over trajectories with period around T, and D is a constant. 
For the evaluation of (54) we have to evaluate the averages 



j(fl7-y) 



(56) 

for the different kinds of correlated orbit pairs. In the diagonal approximation one pairs orbits 
that are either identical or traverse the same path in the opposite direction. In the first case 
the phases are identical and the average in (56) is one. In the second case the phases have 
opposite sign, and with the assumption of the Gaussian distribution of the phases the average 
in (56) becomes 

W=e-* r ^ (57) 



where y = 4Y 2 DTh / h 2 . The remaining evaluation of the diagonal approximation follows section 
4 and the result is 

Te -(B+M)r ( 1+e -yrj (5g) 

It interpolates between the GOE and GUE results as y goes from zero to infinity. 

For the correlated periodic orbit pairs that give the off-diagonal contributions to the form 
factor, we have to distinguish parts of the trajectory that are followed in the same direction and 
parts that are followed in the opposite direction. This depends on the "structure" of the orbit 
pair. Whereas in all previous calculations in this article one could combine the contributions 
from all structures with the same vector v, one now has to look at the structures individually. 
Again one has to treat the orbit stretches and the loops separately. 

Let us look first at the phase differences that arise from the encounter regions. In each 
encounter region a there are l a orbit stretches. Each stretch contributes a value ±0° nc to the 
total phase, and the sign depends on the direction in which the stretch is traversed. The total 
phase difference between an orbit and its partner that comes from an encounter region is then 
±2n a enc where n a is an integer. So for each structure one has to determine the numbers 
n a for all encounters. These numbers can be obtained from the permutation matrices. The 
contribution to the average (56) from an encounter region a is a factor 

^ e ±2in a 0« nc /^ _ e -niyt™ nc /T H 

In order to calculate the phase difference from the loops one has to specify for each structure 
the number of loops iV that are traversed in opposite directions by an orbit and its partner. 
Let ti, . . . , t N denote the times along these N loops. Then the contribution of the loops to the 
average (56) is a factor exp(— y(ti + • • • + t N )/T H ). Of course, the times along the loops vary 
for different orbit pairs which have the same structure. The final result for the average (56) is 

( e ' (fl-y h fly) ) = e -E a <yte n JTH e -y(t 1 +-+t N )/T H ^ 
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In order to evaluate the contribution of orbit pairs with the same structure to the form 
factor, one has to discuss how the density of self-encounters (36) was derived. It was obtained 
by integrating over all possible loop lengths [13] 



w T (s,u) = QL _ vrr - / dt L -! / dt L _ 2 .../ dh (61) 

14 lla r cnc JO JO JO 

where t C nc = J2 a ^enc i s the sum of all the encounter times. The integral is over L — 1 loops, 
the time along the remaining L-th loop is determined by the condition that the sum of all loop 
times plus t enc is T. We can now evaluate the effect of the average over the phase differences 
on the contribution to the form factor by including (60) in the integral (61), and we define 

*■<••«>= m-n.it, 

T — tone /"T— t cnc — pT—t ejlc —tL — l--- — t2 

< I dt L _! / dt L _ 2 ... / dtie -y(t 1+ ...+t N )/T H ( 62 ) 



JO Jo 



These integrals can easily be obtained, but they depend now on N as well. 

Finally, the effect of the escape probability in open systems and the variation of the first 
non-symmetry breaking parameter is taken into account in the same way as in the previous 
section. The contribution of a particular structure to the form factor can then be written as 

K stT (r, x,y,M) = rJ d L ~ v s d L ~" 'u Mf^le^ . (63) 

where 

Z T (S, U) = W T (S, U) e ~BT/T He -MT/T H "Q e -(Bl a -M)(l a -l)t ?nc /T H ^ 

a 

The remaining calculation can be done by a computer program. The integrals in (63) are 
treated by using (37). This leads again to the condition that after an expansion of zt(s,u) in 
all encounter times t" nc only those terms survive where all encounter times have a zero exponent. 
All the relevant information about the quantities v, n a and iV for the different structures can 
be found in table 2 of [18]. We list the result in table 3. In order that this table is not too 
complicated we have summed up the contributions to the form factor from all structures with 
the same vector v. We defined y = yr. 

Finally, we sum over all contributions K v (r,x,y, M) in table 3 and add the diagonal ap- 
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Table 3: Contribution of different types of orbit pairs to the form factor K(r,x,y, M). 



proximation. This gives the expansion of the form factor up to order r 5 
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(65) 



This result encompasses all previous results in this article, up to order r 5 . In order to compare 
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it to the result of RMT, we have to expand the exponentials e yr in r as discused in section 3. 
If one does this, one obtains an expansion that agrees with all terms in (22). 

7 Conclusions 

In this article we extended previous work to include open systems. We considered correlation 
functions of the Wigner time delay, which, because it can be expressed in terms of the periodic 
orbits of the repeller, can be treated in a similar way to previous calculations. Compared to 
the closed system case, the calculation must be performed by taking into account the average 
probability of survival of periodic orbits. A small, but important correction comes from the 
encounter regions, where, because the traversals are spatially close, the probability that each 
survives is essentially the probability that the encounter region survives [15]. By incorporating 
this correction we derive the open form factor up to 9th order for the case without time reversal 
symmetry and up to 7th order for the case with time reversal symmetry. These results showed 
exact agreement with the small r expansion of the GUE and GOE results respectively. When 
the number of open channels M is set to zero, the system is closed, and these results agree with 
the terms in the spectral form factor expansion. 

We also took into account the effects of a non-symmetry breaking external parameter. 
The parametric form factor involves comparing the time delay at two diffent values of this 
parameter, correctly rescaled. The semiclassical calculation assumes a Gaussian distribution of 
the parametric velocities of long orbits and again requires a small correction from the encounter 
regions. This correction comes from the fact that the traversals of the encounter region are 
correlated, affecting the Gaussian averaging. We derived the open parametric form factor 
again up to 9th and 7th order for the GUE and GOE cases and found agreement with the 
RMT expansions. When the scaled parameter difference x was set to zero, effectively excluding 
parametric correlations, we recover the results of the form factor for the time delay without an 
external parameter. Likewise, closing the system by setting the number of open channels M to 
zero allows us to recover previous parametric results [18, 19]. 

The final, and most general, case we considered also allowed for a second symmetry breaking 
parameter y, describing a magnetic field. This parameter gives an additional term that depends 
on the time an orbit and its partner spend travelling in opposite directions. Semiclassically, 
the calculation is complicated by the need to consider each "structure" separately, as it now 
depends on the number of loops and encounter traversals that are followed in opposing directions 
[16, 18]. We obtain a semiclassical result up to 5th order in r, from which we can recover the 
parametric form factor for the time delay by taking the limit y = to get the GOE case and 
y — > oo to get the GUE case. For fixed y we recover the small r expansion of the RMT result. 
Also by removing the parametric correlation x — 0, and closing the system M = 0, we recover 
the GUE-GOE transition result [16, 18]. 

The semiclassical calculation for all these cases agrees with the small r expansion of the 
appropriate RMT result. The limitations of the calculation are similar to that for the spectral 
form factor. A main difficulty is to show that there are no other semiclassical contributions that 
survive the semiclassical limit and would lead to a deviation from the RMT results. We also 
assume that the average over survival times and Gaussian distributed parametric velocities can 
be done indepentently of the sum over orbits. Another open point concerns the region r > 1, 
in which the random matrix expressions for K(r, x, y, M) have a different functional form. 
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